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This  paper  reports  on  the  design,  fabrication  and  proof  of  concept  of  a  multilayer  fluidic  packaging  sys¬ 
tem  enabling  an  increase  in  the  output  power  performance  of  micro  thermoelectric  generators  (pTEGs). 
The  complete  integration  of  the  microfluidic  heat  transfer  system  (fjiHTS)  with  a  pTEG  is  successfully 
demonstrated. 

The  fabricated  prototype  is  characterized  with  respect  to  its  thermal  and  hydrodynamic  performance 
as  well  as  the  generated  output  power.  At  a  very  low  pumping  power  of  0.073  mW/cm2,  a  heat  transfer 
resistance  of  0.74  cm2  K/W  is  reached.  The  assembled  device  generated  up  to  1.47  mW/cm2  at  an  applied 
temperature  difference  of  50  K  and  a  fluid  flow  rate  of  0.1 1/min.  Further  system  improvements  and  the 
potential  of  the  proposed  packaging  approach  are  discussed. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

With  significant  improvements  in  the  figure  of  merit  (ZT)  of 
thermoelectric  materials  [1,2]  and  the  potential  of  large  scale 
production  [3],  thermoelectric  energy  harvesting  is  becoming 
promising  for  commercial  applications.  Recovered  waste  heat  from 
industrial  equipment,  automotives  or  home  appliances  can  be  used 
to  power  a  variety  of  systems,  such  as  wireless  sensor  nodes  [4],  car 
electronics  [5]  or  small  consumer  applications. 

Efficient  thermoelectric  power  generation,  however,  requires 
not  only  a  high  performance  TEG  but  also  a  complete  system 
optimization  including  thermal  coupling  to  the  hot  and  cold  side 
reservoir.  The  output  power  of  a  generator  scales  with  the  square  of 
the  applied  temperature  gradient  across  the  TEG.  In  order  to  main¬ 
tain  this  thermal  gradient,  efficient  heat  dissipation  and  supply  are 
essential,  in  particular  for  thin  pTEGs. 

The  most  common  approaches  to  maintain  a  thermal  gradient 
across  a  TEG  are  the  integration  of  passive  heat  spreaders  or  con¬ 
ventional  large  scale  fluidic  heat  exchangers  [6].  The  first  solution 
has  the  advantage  of  not  requiring  an  active  cooling  medium,  but 
suffers  from  high  heat  transfer  resistance  and  thus  insufficient 
heat  transfer.  Large  scale  fluidic  heat  exchangers  have  an  improved 
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thermal  performance,  but  also  consume  much  pumping  power 
due  to  high  fluid  flows. 

A  more  effective  approach  to  enhance  the  thermal  coupling  to 
the  hot  and  cold  reservoirs  is  to  increase  the  convective  surface 
area  by  using  heat  transfer  systems  on  the  micrometer  scale  [7]. 
Another  approach  is  to  reduce  the  interfacial  thermal  losses  by 
direct  integration  of  heat  transfer  systems  (HTSs)  into  the  pack- 
aging  [8], 

Many  efficient  solutions  for  active  surface  cooling  in  the 
micrometer  range  have  been  developed  for  integrated  circuits 
and  optoelectronic  devices.  The  system  designs  and  methods 
include  simple  microchannels  [9],  multilayer  manifold  structures 
[10-12],  porous  media  [13]  and  jet  impingement  [14].  While  most 
devices  exhibit  a  similar  heat  dissipation  capacity,  manifold  struc¬ 
tures  in  particular  show  a  superior  hydrodynamic  performance 
[7]. 

A  representative  comparison  of  the  impact  of  different  heat 
dissipation  systems  on  the  output  performance  of  a  Bi2Te3  pTEG 
(Z=  1.1  x  10-3  K-1)  is  depicted  in  Fig.  1.  The  heat  dissipation  sys¬ 
tem  performance  is  characterized  by  its  heat  transfer  resistance 
Rhs-  The  pTEG  output  power  is  simulated  as  a  function  of  differ¬ 
ent  thermal  coupling  to  the  cold  side  reservoir  (T0)  and  a  varying 
heat  source  temperature  (Thot)-  The  applied  model  consists  of  a 
TEG  model  with  matched  electrical  load  [15]  extended  with  an 
attached  heat  transfer  resistance,  Rhs-  Fig.  1  shows  that  a  signif¬ 
icant  increase  in  the  output  power  can  be  achieved  by  reducing 
the  heat  transfer  resistance.  Therefore,  combining  pTEGs  with  low 
pumping  power  microfluidic  HTS  will  result  in  a  substantial  output 
power  enhancement. 
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Fig.l.  Simulated  Bi2Te3  |jlTEG(Z=1.1  x  10-3  K_1,J?Teg  =  5.9  cm2  K/W),  output  power 
dependence  on  different  attached  heat  transfer  systems  on  the  cold  side  with  repre¬ 
sentative  values  of  heat  transfer  resistances  (J?hs)-  The  cold  side  temperature  is  fixed 
to  To  =  300  K.  A  perfect  heat  transfer  between  7Hot  and  the  pTEG  was  assumed.  The 
figure  inset  shows  the  schematics  of  the  simulated  setup  and  its  equivalent  thermal 
circuit. 


two-layer  configuration  results  in  very  short  fluid  paths  inside  the 
microchannels,  ensuring  uniform  cooling  and  small  pressure  losses 
[10].  In  particular,  the  reduced  hydrodynamic  resistance  offers  a 
significant  advantage  over  conventional  one  layer  heat  sinks,  since 
the  necessary  pumping  power  needs  to  be  subtracted  from  the 
pTEG  output  power  for  an  overall  system  performance  evaluation. 
Additionally,  low  heat  transfer  resistances  and  thus  high  heat  fluxes 
are  enabled  by  the  high  aspect-ratio  copper  microchannels  and 
the  thin  interface  bond  layer  between  [xHTS  and  pTEG  (Fig.  2a).  A 
further  system  improvement  is  achieved  by  choosing  a  thermally 
isolating  polymer  as  the  manifold  material.  This  limits  the  ther¬ 
mal  crosstalk  between  the  manifold  and  the  copper  channels  and 
therefore  the  heating  up  of  the  inlet  fluid. 

In  this  work,  the  fabrication  process  of  two  layer  fxHTSs  and 
their  integration  with  pTEGs  is  presented.  A  fabricated  p,HTS  is 
first  analyzed  with  respect  to  its  heat  transfer  and  hydrodynamic 
performance.  In  a  second  step,  a  |xHTS  device  is  integrated  with 
a  pTEG  and  the  thermoelectric  Seebeck  voltage  is  measured  as 
a  function  of  the  applied  fluid  flow  rate  and  temperature  gradi¬ 
ent.  Based  on  the  measurement  results,  the  influence  of  the  |xHTS 
to  pTEG  thermal  resistance  ratio  as  well  as  the  potential  per¬ 
formance  improvement  of  the  microfluidic  package  approach  is 
discussed. 


2.  Fabrication  process 


Potential  application  fields  of  pTEGs  with  fluidic  heat  dissipa¬ 
tion  and  supply  would  range  from  automotive  systems  to  industrial 
processes,  where  cold  as  well  as  warm  fluids  are  available. 

The  presented  packaging  approach  enables  the  integration  of  a 
multilayer  fluidic  |xHTS  by  directly  bonding  it  onto  the  pTEG  sur¬ 
face.  In  Fig.  2a,  an  exploded  view  of  the  pTEG/p.HTS  stack  is  depicted 
and  Fig.  2b  shows  the  assembled  device.  The  packaging  system  con¬ 
sists  of  a  high  aspect-ratio  copper  microchannel  layer  as  the  heat 
transfer  structure  and  a  polymeric  manifold  layer.  The  manifold 
layer  homogenously  distributes  the  liquid  through  the  inlet  slots 
into  the  lower  microchannels  and  collects  it  at  the  outlet  slots,  as 
schematically  shown  in  Fig.  2c.  For  a  constant  pressure  drop  over 
the  entire  channel  length,  the  manifold  channels  are  tapered.  This 
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Fig.  2.  Exploded  view  of  the  (jlHTS  and  pTEG  multilayer  stack  (a).  Assembled  system 
(b)  and  cross  section  view  showing  the  fluid  path  (c). 


The  two  functional  layers  of  the  p-HTS,  the  microchannel 
and  manifold  structure,  respectively,  are  fabricated  in  separate 
processes  and  assembled  on  chip  level.  For  the  microchannel  fab¬ 
rication,  an  SU-8  based  LIGA  (Lithography,  Electroplating,  and 
Molding)  process  was  used,  allowing  for  high  aspect  ratios  [16,17], 
batch  fabrication  and  low  cost  electrochemical  metal  deposition. 
The  process  flow  is  illustrated  in  Fig.  3.  In  the  first  step,  1  pm  of  lift¬ 
off  resist  (LOR)  is  spin  coated  on  the  wafer  for  a  later  release  of  the 
microchannels  from  the  support.  Next,  a  copper  seed  layer  (1  pm) 
and  a  chromium  adhesion  layer  (10  nm)  are  deposited  by  e-beam 
evaporation.  The  prepared  seed  wafer  is  coated  with  a  200  pm  thick 
SU-8  3050  film.  To  improve  the  layer  thickness  homogeneity,  the 
spin  coating  is  performed  in  two  steps  of  100  pm  and  an  edge  bead 
removal  is  carried  out.  In  the  subsequent  photolithography  (Fig.  3a), 
long  cooling  and  additional  relaxation  times  are  added  for  stress 
reduction  in  the  SU-8.  Due  to  the  high  aspect  ratio  and  length  of  the 
mold  structures,  a  reduced  intrinsic  stress  level  of  SU-8  is  crucial  for 
the  adhesion  to  the  substrate.  The  structures  are  developed  in  ultra¬ 
sound  in  order  to  improve  convection  inside  the  high  aspect  ratio 
trenches.  The  copper  microchannels  are  grown  by  electrochemical 
deposition  with  an  applied  current  density  of  3  A/dm2  (Fig.  3b).  The 
mold  is  overgrown  in  order  to  form  the  base  of  the  microchannels. 
After  a  top  surface  planarization  by  mechanical  polishing,  the  wafer 
is  released  from  the  support  and  separated  into  8  mm  x  8  mm  chips. 
The  seed  layer  peels  off  during  the  release  process.  Finally,  the  SU-8 
mold  is  removed  by  chemical  etching  (Fig.  3c). 

The  microchannels  are  bonded  to  the  pTEG  (the  basic  pTEG  fab¬ 
rication  process  is  reported  in  [3])  by  an  8  pm  benzocyclobutene 
(BCB)  adhesion  layer  (Fig.  3d).  The  bonding  is  performed  at  a  low 
temperature  of  423 1<  for  16  h  and  a  tool  pressure  of  30kPa  in 
order  not  to  damage  the  pTEG.  The  PMMA  manifold  distribution 
channels  are  fabricated  by  micromilling  and  encapsulated  by  7- 
butyrolactone  (GBL)  solvent  bonding  of  a  250  pm  PMMA  foil.  A 
10  pm  thin  double-sided  adhesive  tape  is  laminated  to  the  man¬ 
ifold,  and  the  connecting  slots  at  the  bottom  of  the  manifold 
structure  (see  Fig.  2c)  are  opened  by  micro  laser  milling  (Fig.  3e). 
Finally,  the  manifold  is  mounted  onto  the  microchannel  and  the 
fxTEG  stack  (Fig.  3f). 
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3.  Experimental 

The  fabricated  |jlHTS  and  p,HTS/p,TEG  stack  are  characterized 
by  heat  flux,  pressure  drop  and  Seebeck  voltage  measurements  in 
a  fluid  loop  measurement  setup.  The  devices  are  mounted  on  a 
temperature-controlled  heat  source,  consisting  of  a  resistive  heater 
(platinum  meander  on  A1203  substrate)  and  a  BCB-bonded  copper 
block  for  homogenous  heat  distribution  (see  inset  in  Figs.  5  and  6). 
The  Cu  block  temperature  is  controlled  with  the  help  of  an  inte¬ 
grated  thermocouple  located  at  the  center  of  the  block.  The  surface 
area  of  the  Cu  block  exposed  to  ambient  air  is  kept  small  in  order 
to  reduce  thermal  losses  through  convection.  Thus,  the  convective 
thermal  losses  are  neglected.  For  an  improved  thermal  interface 
between  the  Cu  block  and  the  heat  sinks,  a  silicone  thermal  com¬ 
pound  (k  =  3.4W/mI<)  is  used  and  a  constant  pressure  of  230  kPa 
applied.  The  heat  flux  through  the  device  is  assumed  to  be  equiva¬ 
lent  to  the  necessary  heating  power  of  the  resistive  heater  in  steady 
state  condition.  The  cooling  water  (deionized)  flow  is  controlled 
by  a  Coriolis  flow  meter  (error:  ±1%  of  reading)  from  0.015  to 
0.13 1/min.  The  water  inlet  temperature  is  kept  at  293  ±  0.1  K.  The 
applied  temperature  difference  AT  is  calculated  by  subtracting  the 
water  inlet  temperature  from  the  controlled  heater  block  temper¬ 
ature.  The  pressure  drop  over  the  heat  sink  is  measured  with  a 
differential  pressure  sensor  meter  (error:  ±2%  of  FS  =  1  mbar).  The 


Fig.  4.  Pictures  of  the  assembled  |jlHTS/(jlTEG  stack  (a)  and  cross-section  view  of  a 
representative  sample  (b).  The  microchannels  are  filled  with  copper  residues  from 
the  polishing  process  for  cross-section  preparation.  Top-view  of  a  microchannel  chip 
and  the  corresponding  cross  section  (c). 

heat  flux  and  pressure  loss  measurements  are  sampled  at  a  rate  of 
1.17  Hz  and  averaged  over  3  min  after  reaching  steady  state  condi¬ 
tion.  The  |jiHTS  performance  is  characterized  by  the  achieved  heat 
transfer  resistance  (RHs)  and  the  corresponding  pumping  power. 
The  heat  transfer  resistance  is  calculated  from  the  measured  heat 
flux  and  the  applied  temperature  difference  AT.  The  pumping 
power  is  obtained  from  the  resulting  pressure  drop  over  the  heat 
sink  and  the  applied  volumetric  flow  rate. 

For  the  p^HTS/pTEG  system  characterization,  the  open  circuit 
voltage  of  the  p/TEG  is  recorded  as  a  function  of  the  fluid  flow  rate 
and  the  applied  temperature  gradient.  The  relevant  parameters  of 
the  |jlTEG  used  can  be  found  in  Table  1. 

4.  Results  and  discussion 

4.  t .  Device  fabrication 

An  assembled  ixHTS/pTEG  system  is  depicted  in  Fig.  4a.  The 
cooling  liquid  enters  the  manifold  on  the  left  hand  side,  flows 
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Table  1 

Bi2Te3  pTEG  parameters. 


Seebeck  voltage  per  TC 
(pA//K) 

Efficiency  ZT(-) 

TEG  output  power 
(|jlW/K2  cm2) 

TEG  resistance  (£2) 

Total  thickness  (pm) 

Thermal  TEG  resistance 
(cm2  K/W) 

102.3  ±0.76 

0.01 64  ±0.001 

0.95  ±0.0139 

5.2  ±6.12E-4 

290  ±7.1 

9.13  ±0.41 

Table  2 

Microchannels  and  manifold  design  parameters. 


Channel  height 
(pan) 

Channel  width 
(pm) 

Fin/channel  wall 
width  (pm) 

Base  height  (pm) 

Slot  size  (pm) 

Total  area  (mm) 

Cu  microchannels 

200 

35 

35 

10 

- 

00 

X 

00 

Manifold  channels 

1000 

200-390  (tapered) 

120 

200 

150 

8.5x12 

Fluid  Flow  [l/min] 


Fig.  5.  Heat  transfer  resistance  Rhs  and  corresponding  pumping  power  measure¬ 
ments  and  calculations  of  the  pHTS  as  a  function  of  the  fluid  flow  rate.  The  Cu  block 
temperature  was  controlled  to  298  ±0.1 1<.  The  error  bars  correspond  to  the  stan¬ 
dard  deviations  of  the  averaged  signal  over  a  measurement  period  of  180  s.  The 
figure  inset  shows  the  schematics  of  the  sample  fixture  in  the  measurement  setup. 


through  the  microchannels  located  underneath  and  exits  the  man¬ 
ifold  on  the  opposite  side.  The  device  is  connected  to  the  fluid 
loop  by  a  flexible  tube  at  the  fluid  inlet  and  outlet.  The  relevant 
device  parameters  are  summarized  in  Table  2.  The  dimensions  of 
the  |jiHTS  were  chosen  to  maximize  the  output  power  according 
to  an  extended  thermal  resistance  model  of  the  pHTS/pTEG  stack 
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Fig.  6.  Seebeck  voltage  measurements  of  the  pHTS/pTEG  stack  as  a  function  of  the 
fluid  flow  rate  at  an  applied  AT  of  30  ±  0.1  K.  The  error  coming  from  random  voltage 
fluctuations  is  below  0.1  mV/cm2  (not  included  in  graphic).  The  figure  inset  shows 
the  schematics  of  the  sample  fixture  in  the  measurement  setup. 
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[10,15]  while  still  complying  with  fabrication  limits.  Increasing  the 
height  of  both  channels  would  result  in  an  improved  hydrodynamic 
performance  and  reduced  heat  transfer  resistance  in  the  case  of 
the  microchannels,  but  would  also  complicate  the  fabrication  pro¬ 
cess.  Fig.  4b  shows  the  cross  section  of  the  assembled  stack  with  a 
Bi2Te3  ijlTEG  pile,  copper  microchannel  and  a  thin  homogenous  BCB 
bond  layer  in  between.  This  bond  interface  also  serves  as  electric 
isolation  between  the  pTEG  interconnects  and  the  microchannel 
base.  The  depicted  microchannels  and  the  pTEG  in  Fig.  4b  are  taken 
from  previously  fabricated  samples  and  do  not  correspond  to  the 
actual  measured  device.  Fig.  4c  shows  the  top  view  and  a  cross  cut 
view  of  the  used  microchannels.  The  channels  have  an  aspect  ratio 
of  5.4  at  a  total  height  of  190  pm. 


4.2.  /xHTS  performance 


The  measured  heat  transfer  resistance  (RHs)  and  corresponding 
pumping  power  with  respect  to  the  volumetric  flow  rate  are  plot¬ 
ted  in  Fig.  5.  The  applied  flow  rates  of  0.015-0.13  l/min  are  in  the 
relevant  range  for  the  characterization  of  the  pHTS  under  investi¬ 
gation.  The  Cu  block  temperature  was  maintained  at  298  ±  0.04  K. 
With  rising  flow  rates,  RH s  decreases  due  to  an  increase  in  the  heat 
capacity  rate  of  the  fluid.  The  total  thermal  resistance  of  the  pHTS 
is  defined  as 


Rhs  =  Rc  ond  +  ^conv  +  ^hc 


(1) 


where  Rc ond  is  the  conductive  resistance  from  the  heater  to  the 
copper  channel  surface  area,  RConv  the  convective  resistance  and 
Rue  the  fluid  resistance  due  to  its  limited  heat  capacity.  RHc  can  be 
expressed  as 


1 

vpep 


(2) 


where  m  is  the  mass  flow  rate,  V  the  volumetric  flow  rate,  p  the 
fluid  density  and  cp  the  specific  heat  capacity.  While  the  conductive 
and  convective  heat  transfer  resistances  are  nearly  independent  of 
the  flow  rate,  Rhs  decreases  with  increasing  flow  rates  due  to  an 
increase  in  the  effective  thermal  mass  flow.  At  higher  flow  rates,  the 
impact  of  the  enhanced  mass  flow  decreases  and  a  further  rise  in  the 
flow  rate  only  marginally  changes  the  total  heat  transfer  resistance. 

Rc ond  includes  the  thermal  resistance  of  the  copper  channels  as 
well  as  the  interface  losses  between  the  Cu  block  and  the  pHTS. 
While  the  first  can  be  neglected  due  to  the  high  thermal  conductiv¬ 
ity  of  copper,  the  latter  might  include  the  effects  of  a  thick  silicone 
thermal  compound  or  interface  voids,  resulting  in  a  significant 
impact  on  the  overall  heat  transfer  resistance. 

The  pumping  power  in  Fig.  5  increases  with  increasing  fluid  flow 
due  to  a  rise  in  the  hydrodynamic  resistance  with  higher  water  flow. 
In  the  applied  flow  rate  range  (0.015-0.13  l/min),  the  calculated 
Reynolds  number  at  the  manifold  channel  inlet  varies  from  40  to 
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370  and  2  to  20  inside  the  microchannels,  indicating  a  laminar  flow 
regime. 

The  applied  model  in  Fig.  5  is  a  semi-empirical  one-dimensional 
model  [18]  extended  with  temperature  dependent  fluid  parame¬ 
ters.  The  clear  downward  offset  of  the  modeled  ftHs  with  respect 
to  the  measurements  indicate  the  above  mentioned  thermal  losses 
between  the  Cu  block  and  the  fjiHTS  (not  included  in  the  applied 
model).  These  thermal  losses  correspond  to  a  silicone  thermal  com¬ 
pound  thickness  of  approx.  26  |jim.  The  computed  pumping  power 
is  in  good  agreement  with  the  experimental  results. 

Since  the  consumed  power  must  be  subtracted  from  the  |xTEG 
output  power,  a  trade-off  between  the  heat  transfer  resistance 
and  the  pumping  power  must  be  made.  Depending  on  the  spe¬ 
cific  applied  temperature  gradient  and  the  [xTEG  parameters,  an 
optimal  fluid  flow  operating  point  will  exist  (will  be  discussed  in 
Section  4.3). 


4.3.  fiTEG//nHTS  performance 


The  open  circuit  Seebeck  voltage  measurements  as  a  function 
of  the  fluid  flow  rate  through  the  |xHTS  are  depicted  in  Fig.  6. 
With  increasing  fluid  flow  rates,  the  Seebeck  voltage  increases, 
following  the  inverse  trend  of  the  heat  transfer  resistance  behav¬ 
ior  (Fig.  5).  However,  it  can  be  noticed  that  a  rising  fluid  flow 
rate,  i.e.  decreasing  RH s,  has  a  relatively  small  impact  on  the 
Seebeck  voltage  compared  to  the  absolute  measured  value.  The 
reason  lies  in  the  large  difference  between  the  thermal  resis¬ 
tance  of  the  ijiTEG  (#rEG  =  9.13  cm2  K/W)  and  the  ixHTS  range 
(Khs  =  0-1 8-0.73  cm2  K/W). 

The  Seebeck  voltage  of  a  TEG  can  be  expressed  as  [19] 

V  =  maATjEQ  =  mot- - - TAMB  (3) 

^TEG  +  ^Contact 

where  m  is  the  number  of  the  thermocouples,  a  the  Seebeck 
coefficient  of  the  thermoelectric  material,  A Tteg  the  temperature 
gradient  at  the  generator  surface,  ^contact  the  sum  of  the  thermal 
contact  resistances  on  the  hot  and  cold  sides  and  ATamb  the  ambi¬ 
ent  thermal  gradient.  In  the  system  investigated  here,  ATamb  is 
equal  to  the  applied  temperature  gradient  AT. 

By  introducing  the  ratio  X  between  ^contact  and  RjEG  and  assum¬ 
ing  mot  ATamb  to  be  a  constant  which  is  denoted  by  C,  the  equation 
can  be  simplified  to 


V  =  c _ _ 

Rjeg  +  XRteg 


=  c 


1 


1  +x 


(4) 


With  decreasing  X,  i.e.  decreasing  ^contact*  the  Seebeck  voltage 
will  increase.  AtX<l,  i.e.  %eg  > ^Contact,  and  small  variations  of 
^Contact  with  respect  to  ftTEG,  the  impact  on  V  will  be  minor. 

The  same  applies  for  the  output  power  P0ut  of  a  |xTEG.  In  a 
simplified  thermal  resistance  model,  P0ut  can  also  be  described  as 
[20] 


^out  =  C* 


Pteg 


(PtEG  +  ^Contact) 


=  c * 


1 


Pteg(1  +Xf 


(5) 


with 


(6) 


where  T  is  the  mean  temperature  of  the  system,  pc  the  Carnot  and 
Pte  the  thermal  efficiency  of  the  thermoelectric  generator. 

By  reducing  ^contact.  the  output  power  of  the  generator  will 
be  increased  for  arbitrary  PTEG.  Additionally,  for  a  fixed  ^contact* 
the  maximum  of  the  equation  is  reached,  when  PTEG  is  equal  to 
^Contact-  Therefore,  by  decreasing  PTEG  (by  e.g.  reducing  the  gen¬ 
erator’s  thickness)  to  match  the  ^contact,  the  impact  of  the  small 
PHs  can  be  better  exploited  and  the  output  power  can  be  further 


Fig.  7.  Seebeck  voltage  measurements  and  maximum  output  power  of  the 
(xHTS/pTEG  stack  as  a  function  of  applied  temperature  gradient  at  a  fluid  flow  rate  of 
0.1 1/min.  The  error  coming  from  random  voltage  fluctuations  is  below  0.1  mV/cm2 
(not  included  in  graphic). 


enhanced.  However,  in  real  applications,  practical  limits  of  main¬ 
taining  the  temperature  gradient  ATamb  with  small  PTeg  exist  and 
must  be  considered  for  the  system  design. 

In  Fig.  7,  the  measured  Seebeck  voltage  and  the  correspond¬ 
ing  calculated  |xTEG  output  power  is  plotted  versus  the  applied 
temperature  gradient.  A  fluid  flow  rate  of  0.1 1/min  was  selected 
to  explain  the  dependence  of  the  Seebeck  voltage  on  AT  for  a 
characteristic  working  point.  The  solid  line  represents  the  See¬ 
beck  voltage  resulting  from  the  previously  characterized  Seebeck 
coefficient  of  the  thermocouples  (see  Table  1).  At  a  specific  See¬ 
beck  voltage,  the  actual  temperature  gradient  at  the  surface  of  the 
generator  can  be  determined  (see  vertical  dashed  line  in  Fig.  7). 
Assuming  that  an  equal  heat  flux  runs  through  the  generator  at  the 
same  measured  Seebeck  voltage,  the  total  thermal  contact  resis¬ 
tance  between  the  reservoirs  of  the  cold  and  hot  sides  is  calculated 
to  be  2.27  ±0.13  cm2  K/W  (at  AT=30K).  By  subtracting  the  cor¬ 
responding  Rhs  of  0.19  ±  IE-3  cm2  K/W  (see  Fig.  5  at  a  fluid  flow 
rate  of  0.1 1/min),  the  thermal  interface  losses  are  estimated  to  be 
2.08  ±  0.13  cm2  K/W.  This  relatively  high  value  mainly  results  from 
the  insufficient  thermal  coupling  between  the  generator  and  the 
heater  block,  since  the  p,TEG  surface  is  not  totally  flat  and  contains 
some  defects  and  voids. 

The  maximum  output  power  is  computed  from  the  measured 
Seebeck  voltage  and  the  electrical  resistance  of  the  |xTEG  (see 
Table  1)  assuming  matched  load  conditions  [15,21].  A  maximum 
power  of  1.47mW/cm2  is  reached  at  an  applied  temperature  dif¬ 
ference  of  50  K.  At  this  high  flow  rate  of  0.1 1/min  and  thus  high 
pumping  power  of  3.74  mW/cm2,  no  positive  net  output  power  can 
be  achieved  with  the  low  ZT  values  of  the  fxTEG  used. 

In  order  to  predict  the  potential  performance  of  the  presented 
p,HTS/|jiTEG  system,  the  measured  |xHTS  parameters  and  improved 
IjiTEG  thermoelectric  material  properties  are  inserted  into  the 
model  of  Fig.  1. 

The  [xHTS  parameters  inserted  into  the  model  are  selected 
in  order  to  maximize  the  computed  net  output  power  for  an 
applied  temperature  gradient  of  0-65  K.  A  heat  transfer  resis¬ 
tance  of  0.74±6E-3  cm2  K/W  at  a  fluid  flow  rate  of  0.01 5 1/min 
is  used  and  the  corresponding  measured  pumping  power  of 
0.073  ±  lE-3mW/cm2  is  subtracted  from  the  |xTEG  output  power. 
The  estimated  thermal  interface  losses  between  heater  and  device 
are  also  included  in  the  model.  The  simulated  output  power  as  a 
function  of  the  applied  temperature  gradient  and  different  ZT  val¬ 
ues  is  plotted  in  Fig.  8. 
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Fig.  8.  Simulated  Bi2Te3  (jlTEG  (RTEg  =  9.13  cm2  K/W)  output  power  dependence  on 
the  applied  temperature  gradient  and  different  ZT  values.  A  heat  transfer  resistance 
of  0.74  cm2  K/W  and  a  pumping  power  of  0.073  mW/cm2  were  used,  corresponding 
to  a  |jiHTS  fluid  flow  of  0.015 1/min.  The  cold  side  temperature  is  set  to  293  K. 

The  optimal  operating  point  for  this  p,HTS/|jiTEG  configuration 
and  boundary  conditions  is  at  a  low  flow  rate,  since  the  increase  of 
the  required  pumping  power  overweighs  the  limited  impact  of  the 
reduced  s  at  higher  flow  rates.  By  decreasing  ftTE G  (i.e.  increasing 
the  ratio  of  fins  to  Rj eg)  and  the  interface  contact  resistance  to  the 
heater,  the  influence  of  a  small  s  will  increase  and  the  optimal 
operating  point  will  shift  to  larger  fluid  flows. 

The  continuous  red  line  in  Fig.  8  represents  the  power  output  of 
the  |jiHTS/|jiTEG  stack  for  the  ZT  of  the  generator  used  (see  Table  1 ). 
A  positive  net  power  can  be  reached  at  a  temperature  gradient  of 
~12K.  By  increasing  the  ZT  to  a  state-of-the-art  value  of  0.5  for 
Bi2Te3,  a  positive  output  power  starting  from  an  applied  AT  of 
~2.5  K  could  be  achieved. 

5.  Conclusions 

A  new  microfluidic  packaging  approach  for  the  power  enhance¬ 
ment  of  fjiTEGs  is  presented.  The  two  level  heat  transfer  system 
consists  of  high  aspect  ratio  copper  microchannels  for  efficient  heat 
dissipation  and  a  polymer  manifold  fluid  distribution  layer  for  a 
reduced  hydrodynamic  resistance. 

The  fabrication  process  of  the  p^HTS  and  the  integration  with 
a  fxTEG  were  successfully  demonstrated.  The  fabricated  prototype 
was  characterized  by  heat  flux,  pressure  losses  and  Seebeck  volt¬ 
age  measurement  in  a  controlled  fluid  loop.  The  |xHTS  achieved 
a  heat  transfer  resistance  of  0.74  cm2  K/W  at  a  very  low  pumping 
power  of  0.073  mW/cm2.  Further,  a  maximum  [jiTEG  output  power 
of  1.47  mW/cm2  was  reached  at  an  applied  temperature  difference 
of  50  K  and  a  fluid  flow  rate  of  0.1 1/min  (3.74  mW/cm2  pumping 
power). 

The  complete  system  can  be  further  optimized,  by  reducing  the 
pTEG  thickness  and  interfacial  thermal  losses.  Simulations  based 
on  the  measured  system  parameters  but  with  a  state-of-the-art  ZT 
value  of  0.5  for  Bi2Te3,  show  a  net  output  power  gain  starting  from 
an  applied  temperature  gradient  as  low  as  2.5  K. 
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